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Abstract:  A  detailed  static  atomistic  model  of  the  dense,  glassy  polycarbonate 
of  4,4'-isopropylidenediphenol  ('bisphenol-A  polycarbonate',  PC)  is  simulated 
using  a  well  established  technique  that  previously  proved  successful  for 
simple  vinyl  polymers.  Initial  chain  conformations,  which  are  generated 
using  a  Monte  Carlo  technique  including  periodic  continuation  conditions,  are 
relaxed'  using  a  potential  energy  minimization.  Two  sizes  of  microstructures 
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at  densities  of  i.2TJ  g/cm^  were  obtained,  one  with  cube  edge  length  of  18  A  and 
the  other  with  an  edge  length  of  30  A.  Detailed  analysis  of  the  minimized 
structures  indicate  intermolecular  packing  influences  create  a  large  variation 
of  chain  conformations  different  from  the  purely  intramolecular  ground 
states,  but  that  it  is  the  intramolecular  energy  contributions  that  determine 
which  combinations  can  or  cannot  occur.  The  systems  are  amorphous, 
exhibiting  random  coil  behavior,  with  some  weak  intermolecular  correlations 
that  exist  on  a  very  small  scale. 
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I.  INTRODUCTION 

There  has  been  considerable  interest  in  the  molecular  level  structure  of 
the  glassy  polycarbonate  of  4,4'-isopropylidenediphenol  (’bisphenol-A 
polycarbonate',  PC)  and  its  link  to  bulk  mechanical  and  transport  properties. 
Several  investigations  have  been  done  using  both  experimental1  * 1 5 
theoretical^*  methods  to  probe  the  behavior  of  glassy  PC  at  the 

atomistic  scale,  and  whereas  the  experimental  studies  have  been  extensive  and 
precise,  the  theoretical  modelling  to  date  has  not  been  as  rigorous,  and  a 
detailed  molecular  level  model  of  glassy  PC  was  lacking. 

Models  of  amorphous  polymers  at  glass  densities  with  chain  backbones 
less  complicated  than  PC,  mainly  polyethylene  and  other  vinyl  polymen.  have 
been  devised  using  both  static *7. 18  and  deterministic  molecular  dynamic 
(MD)1^'^1  methods  and  have  resulted  in  varying  degrees  of  success.  The 
extension  of  these  techniques  to  more  complex  polymers  has  been  slow.  Models 
of  glassy  polymers  with  chain  backbones  as  complex  as  PC  have  frequently 
resorted  to  gross  simplifications  of  the  structure,  limiting  their  application.  It 
was  our  intent  to  use  a  well  established  simulation  method  that  has  been 
successful  for  simple  vinyl  polymers  and  extend  it  to  PC.  We  present  here  a 
model  of  glassy  PC  using  the  molecular-static  method  first  developed  by 
Theodorou  and  Staler  for  atactic  polypropylene17  (PP)  and  then  modified  for 
atactic  polyvinylchloride  (PVC)  by  Ludovice  and  Suter1^.  It  is  presently  the 
most  detailed  model  of  PC  available  and  has  proven  very  useful  in  the 
accompanying  atomistic  level  study  of  conformational  alterations  of  molecular 
segments  of  potential  importance  in  mechsnical  relaxations  in  PC.^2  we  «hal! 
refer  to  this  study  below  as  III. 
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II.  PREVIOUS  WORK 


The  molecular  level  structure  of  PC  has  been  well  characterized  using 
various  experimental  techniques,  among  which  are  included:  light 
scattering,1  x-ray  diffraction,^'^  and  neutron  scattering7*^  as  well  as  Nuclear 
Magnetic  Resonance  (NMR).1®*1^  A  brief  overview  of  these  studies  follows. 

Light  scattering,1  small  angle  x-ray  scattering,^  and  small  angle 
neutron  scattering7-®  studies  have  shown  that  the  single  chain  behavior  of  PC 
in  the  bulk  glass  demonstrates  random  coil  behavior.  In  addition,  the  wide 
angle  light  scattering  analysis  of  Dettenmaier  and  Kausch1  indicate  single 
phase  behavior  with  no  large  regions  of  strong  intermolecular  orientation 
correlations.  Wide  angle  x-ray^*^  and  neutron  scattering^  have  given  more 
detailed  insight  into  the  atomistic  scale  structure  of  PC.  The  results  of  these 
studies  have  generally  been  interpreted  as  indicating  that  numerous  very 
small  regions  of  enhanced  order1  exist  in  the  bulk,  consisting  of  no  more  than 
2  or  3  chemical  repeat  units.  These  regions  are  believed  to  consist  of  axially 
aligned  chain  segments  in  which  the  carbonate  group  assumes  a  trans.trans 
conformation.  The  regions  are  too  small  to  alter  the  overall  random  coil 
behavior  of  the  chain  or  to  show  large  fluctuations  in  the  optical  anisotropy 
of  PC. 

Various  solid  state  NMR  investigations  have  been  used  to  probe  the 
atomic  level  structure  by  studying  the  chain  mobility.  The  characteristic 

chain  motions  in  glassy  PC  have  very  broad  frequency  distributions, 
indicating  large  site  variations  due  to  differences  in  intermolecular  packing, 
rienrichs  and  Nicely1'  also  have  researched  the  chain  conformations  of  PC  in 
the  amorphous  bulk  and  conclude  that  intramolecular  forcea  largely 


influence  the  preferred  bulk  chain  conformations.  Specifically,  while  the 
phenylene  rings  widely  assume  a  propeller'  conformation  in  the 
isopropylidene  moiety  of  the  repeat  unit,  the  carbonate  group  is  oriented  out 
of  the  plane  of  its  adjacent  phenylene  rings. 

Theoretical  modelling  specifically  on  glassy  PC  has  to  date  been  limited. 
Yannas  and  Luise1^  developed  an  early  model  of  glassy  PC  which  includes  a 
molecular  level  model  that  assumes  hexagonal  packing  of  quasi-lattice  sites 
with  spherical  symmetry.  The  model  includes  intramolecular  interactions 
represented  by  intrinsic  torsion  potentials  and  intermolecular  contribution 
by  Lennard-Jones  type  interactions  between  segments.  Perchak,  J.  Schaefer, 
and  coworkers *2,14  dcvised  a  mode!  consisting  of  a  two  dimensional  lattice  of 
interacting  phenylene  rings  with  varying  degree  of  lattice  flexibility.  Fischer 
et  al.9  in  the  analysis  of  their  neutron  scattering  results  modeled  PC  as 
bundles'  of  parallel  segments,  two  repeat  units  long.  A  detailed  model  of 
glassy,  bulk  PC  has  to  date  been  lacking.  Previous  static17-1®  and  MD^-Zl 
simulations  done  at  glass  densities  have  been  confined  to  studying  polymers 
with  very  simple  backbones.  MD  studies  have  included  models  of  polyethylene 
and  PP,  yet  are  limited  by  available  computer  resources.  However,  static 
models  of  PP  by  Theodorou  and  Suter17  and  of  PVC  Ludovice  and  Suter1®  have 
yielded  excellent  results  while  being  computationally  less  intensive  than  MD 
simulations.  It  is  our  intention  here  to  extend  the  static  simulation  technique 
developed  by  the  aforementioned  authors  to  PC. 
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III.  GENERATION  OF  MODEL  STRUCTURE 


3.1  Model  Parameters 

The  amorphous  PC  microstructures  were  generated  by  the  modeling 
technique  developed  by  Theodorou  and  Suter1?  for  PP  and  extended  to  PVC  by 
Ludovice  and  Suter1**  using  the  force  field  described  in  the  previous 
accompanying  paper  of  this  issue^  to  which  we  will  refer  to  below  as  I.  The 
details  of  the  model  are  as  follows: 

A.  The  model  is  static,  with  temperature  only  entering  with  the  specification 
of  the  density  at  which  the  structure  is  relaxed. 

B.  Spatial  periodic  continuation  conditions  are  employed  for  the  cubical 
simulation  cell. 

C.  The  bond  angles  and  bond  lengths  are  fixed.  Molecular  rearrangement 
occurs  only  through  changing  the  torsion  angles  of  the  bonds  (see  Figure  1). 

D.  Methyl  groups  are  represented  as  pseudo-atoms  for  computational 
expedience. 

E.  A  Lennard-Jones  energy  function  is  used  to  represent  the  nonbonded  van 
der  Waals  interatomic  interactions  while  the  Coulombic  interactions  are 
represented  as  a  Block-Walker  function.  18.24  intrinsic  torsional  functions 
specific  to  the  PC  repeat  unit  are  included.^ ^  Both  the  Lennard-Jones  and  the 
Block-Walker  functions  are  splined  due  to  the  finite  size  of  the  periodic 

cube.1 7*1 8,24 

F.  The  density  is  1.20  g/cm^,  appropriate  of  polycarbonate  at  room 
temperature. 2  5 
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G.  Thirteen  independent  microstructures  have  been  generated  with  periodic 
cube  edge  lengths  of  18.44  A.  degree  of  polymerizations  of  35  (MW  = 

4.532  g/mol),  and  with  485  atoms  or  atom  groups  in  each  cube. 

H.  Two  larger  sized  microstructures  have  also  been  completed,  which  have 
periodic  cube  edge  lengths  of  29.87  A,  degree  of  polymerization  of  151 
(MW  =  19.264  g/mol),  and  contains  2051  atoms  or  atom  groups  in  the  cube. 

3.2  Initial  Guess 

The  generation  of  a  starting  chain  conformation  or  'initial  guess' 
structure  follows  the  'modified  Markov'  Monte  Carlo  process  described  by 
Theodorou  and  Suter^  with  one  major  alteration.  In  the  original  scheme,  the 
chain  is  generated  in  a  sequential  bond-by>bond  fashion,  where  at  each  step 
of  the  generation  short  range  (via  rotational  isomeric  state  (RIS) 
probabilities26)  and  long  range  interaction  energies  are  included  in  the 
weighing  of  the  possible  RIS  conformations  of  that  bond.  For  both  PP  and  PVC, 
the  addition  of  the  long  range  interactions  in  the  conformational  probabilities 
of  the  RIS  state  generated  initial  guesses  that  were  sufficiently  spatially 
uniform  to  be  readily  minimized,  even  though  possible  future  overlaps  of  the 
chain  were  not  anticipated  in  the  generation  of  the  initial  guess  (the  RIS 
probabilities  are  bidirectional,  but  the  inclusion  of  the  long  range  energy  is 
only  unidirectional).  While  this  was  not  a  problem  for  PP  or  PVC.  it  did  prove 
to  be  significant  in  the  generation  of  initial  guess  structures  of  PC.  The 
energies  of  the  initial  guess  structures  generated  using  the  sequential  Monte 
Carlo  scheme  were  on  the  order  of  10®  -  10^  kcal/mol,  with  the  majority  of 
the  structures  in  the  10**  -  101 2  kcal/mol  range.  The  bond  segments  tre  much 
larger  and  bulkier  for  PC  than  for  PP  or  PVC,  and  therefore  the  PC  chain 
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segments  frequently  overlapped  due  to  the  unidirectionality  of  the  generation 
technique.  The  resulting  structures  were  spatially  too  unevenly  distributed, 
and  many  could  not  be  minimized.  Therefore,  the  initial  guess  generation 
routine  was  modified.  Rather  than  considering  the  RIS  states  of  each  bond 
sequentially  in  the  Monte  Carlo  scheme,  a  simultaneous  choice  of  the  RIS 

states  of  .he  three  virtual  bonds  of  the  isopropylidene  group  was  made  (the 

bonds  labelled  1,2.  and  3  in  Figure  1)  and  included  together  as  one  unit,  then 
the  conditional  probabilities  of  each  possible  conformation  of  the  unit  as  a 
whole  were  calculated.  The  structures  generated  with  the  modified  routine  had 
energies  on  the  order  of  10®  -  10* 7  keal/mol  with  the  majority  of  the 
structures  having  energies  in  the  range  of  10®  -  10^  kcal/mol.  More 
importantly,  the  generated  structures  were  more  spatially  uniform  than  the 
structures  generated  using  the  unmodified  routine.  The  root-mean-squared 
end-to-end  distances  (<r7>0/M)  for  100  generated  initial  guesses  were 
calculated,  with  an  average  of  1.0  (±  0.4)  A  mol  g  for  degree  of 

polymerization  X  =  35  and  an  average  of  1.2  (±  0.5)  A  mol  g  for  degree  of 

polymerization  X  =  151.  These  average  values  compare  well  with  the 
experimental  range7*®*77*7®  of  0.75  to  1.28  and  indicate  that  the  technique 
correctly  generates  conformations  exhibiting  random  coil  behavior.  These 
initial  guess  structures  were  then  screened  with  respect  to  their  energy 
(taking  those  of  lowest  energy)  and  <r7>  0/M  for  minimization. 

3.3  Relaxation  to  Mechanical  Equilibrium 

The  energy  minimization  procedure  followed  the  same  principles  as  in 
the  technique  developed  by  Theodorou  and  Suter.*7  The  Eulerian  angles  (the 
angles  that  describe  the  orientation  of  the  chain  with  respect  to  the  frame  of 


reference  of  the  cube)  and  the  bond  torsion  angles  (except  the  first  and  last 
along  the  chain)  are  considered  to  be  the  degrees  of  freedom  for  the  system. 
Analytic  expressions  for  the  total  potential  energy  with  respect  to  the 
Eulertan  and  the  rotation  angles  were  derived.  A  quasi-Newton  matrix- 
updating  algorithm  due  to  Broyden,  Fletcher,  Goldfarb,  and  Shanno^  was  used 
to  minimize  the  structures.  The  initial  guesses  to  be  minimized  were  chosen  in 
regards  to  their  energy  and  so  that  their  average  root-mean-squared  end-to- 
end  distances  agreed  well  with  experimental  values.  The  average  <r2>Q/M  for 

the  smaller  cubes  was  0.84  ±  0.42  (Amolg'1)  and  for  the  larger  cubes  was  0.88 
±  0.09  (Amolg1). 

During  the  minimization,  PC  again  proved  to  be  much  more  demanding 
than  PP  or  PVC.  For  PP,  the  initial  structures  were  minimized  in  three  steps 
using  the  following  modifications  to  the  bonded  and  nonbonded  interactions: 

1) .  The  first  step  included  nonbonded  forces  represented  by  a  soft- 
sphere  function  (which  only  includes  the  repulsive  part  of  the 
Lennard-Jones  interactions),  where  the  atomic  radii  were  taken  at  half 
their  actual  value.  No  rotational  barriers  were  included. 

2) .  In  the  second  step  again  a  soft-sphere  function  was  used,  but  with 
atomic  radii  increased  to  their  actual  size.  Full  rotational  barriers  were 
included. 

3) .  For  the  final  step,  a  full  Lennard-Jones  function  was  introduced  for 
the  nonbonded  interactions  with  atomic  radii  taken  at  their  actual  size. 
Full  rotational  barriers  were  included. 

For  PVC,  with  prominent  Couiombic  charges,  Ludovice  and  Suter  chose  a 
four  step  minimization  which  included  the  following  features:*® 
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1) .  The  first  step  employed  a  soft-sphere  function,  and  the  atomic  radii 
and  Coulombic  charges  were  taken  at  half  their  actual  size.  No 
rotational  barriers  were  included. 

2) .  In  the  second  step,  a  soft-sphere  function  was  used  with  the  atomic 
radii  and  Coulombic  charges  increased  to  70%  of  their  actual  size.  Full 
rotational  barriers  were  included. 

3) .  A  soft-sphere  function  was  again  employed  for  the  third  step,  with 
the  atomic  radii  and  Coulombic  charges  increased  to  their  actual  size, 
and  full  rotational  barriers  included. 

4) .  The  final  step  included  the  full  Lennard-Jones  function,  with  the 
atomic  radii  and  Coulombic  charges  taken  at  their  actual  size,  and  with 
full  rotational  barriers. 

For  PC.  it  was  found  that  many  more  steps  were  necessary,  and  the 
minimization  was  performed  in  two  major  stages: 

Stage  I:  First,  using  only  the  soft  sphere  function  and  full  rotational 
barriers,  the  sizes  of  the  atomic  radii  were  slowly  increased  (stepwise 
by  2.5%  to  10%)  while  simultaneously  the  Coulombic  charges  were 
slowly  increased  to  their  full  values,  for  an  average  of  14  to  20  steps. 

Stage  2:  Then,  with  the  full  Lennard-Jones  functiou  and  full  rotatiooal 
barriers  included,  the  atomic  radii  and  the  Coulombic  charges  were 
again  slowly  increased  to  their  actual  sizes,  over  an  average  of  13  to  20 
steps. 
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Overall,  more  than  25  steps  were  required  in  the  minimization,  which 
is  significantly  more  than  for  e.ther  PP  or  PVC.  The  gradients  o?  energy  for  PC 
are  also  higher  than  foi  PP  or  PVC,  and  for  the  X  =  35  microstructures  were 
less  than  10*3  kcal/fmol  deg)  and  for  the  X  =  151  microstructures  less  than 
10*’  kcal/(mol  deg).  The  difficulty  in  the  minimization  and  the  higher 
energy  gradients  are  most  likely  due  to  the  complexity  of  the  PC  repeat  unit 
and  the  bulkiness  associate  with  the  degrees  of  freedom  in  the  minimization, 

in  Figure  2  a  fully  minimized  structure  fo:  X  =  35  is  shown. 

This  complex  mirimization  procedure  brings  into  question  the  role  of 
the  path  dependence  of  the  minimization  on  the  characteristics  of  the  final 

structures.  To  investigate  this  aspect,  two  different  paths  were  employed.  For  6 
of  the  13  structures  with  degree  of  polymerization  X  =  35  that  were  minimized, 
the  Coulorabic  charges  were  included  only  during  the  second  stage  ot  the 

minimization.  Upon  analysis,  t be  two  sets  of  X  *  35  structures  showed  no 
statistical  differences  form  each  other,  so  all  13  were  considered  as  one 
ensemble. 

IV.  RESULTS 

4,1  Coheaive  Energy  Density 

The  cohesive  energy.  UCoh*  »»  ^e  energy  associated  only  with  the 
intermolecular  forces  and  can  be  estimated  from  an  ensemble  of 
microstructures  by  taking  the  difference  between  the  total  energy  of  the 

microstructure,  Utot>  and  that  of  the  parent  chain  only,  Uptr-  Thus, 


1  1 


Ucoh  *  (Up»r  Utot) 


(D 


which  relates  to  Hildebrand's  solubility  parameter  by  the  following: 

5  =  (UC0h/V)1/2  (2) 

where  V  is  the  volume  of  the  microstructure. 

The  various  energies  were  determined  in  a  similar  manner  to  the 
calculations  preformed  on  PP^7  and  PVC. *8,24  i„  or(jcr  t0  determine  the 
parent  chain  energy  of  the  minimized  structures,  the  cube  edge  length  was 
increased  by  two  orders  of  magnitude  to  fully  enclose  the  chain  in  the  box. 

The  energy  of  this  'isolated'  chain  was  then  taken  as  the  parent  chain  energy. 
Upar- 

In  the  calculation  of  the  total  energy  of  the  microstructure.  Ut0t.  since 
the  nonbonded  energy  interactions  had  been  splined  in  the  minimization,  a 
tail  correction'  factor  had  to  be  introduced.  With  the  inclusion  of  the 
Couiombic  interactions,  the  tail  correction  calculation  becomes  quite  complex, 
leadiag  to  upper  and  lower  bounds  for  the  tail  correction  where  the  lower 
bound  is  achieved  by  assuming  the  structure  is  completely  amorphous  and  the 
upper  bound  is  calculated  assuming  liquid  like  ordering  of  the  dipoles. ^*2 3,2 4 

Overall,  the  cohesive  energy  becomes: 


Ucoh  -  Uper  *  Utot  -  UpM  -  (  U+  AUtlii,  J  (3) 

where  U  is  the  total  energy  of  the  microstructure  and  A  U  tails  **  the  tail 
correction.  Table  1  lists  the  calculated  energies  for  the  thirteen  18  A  cubic 
structures  (degree  of  polymerization  X  ■  33)  and  for  the  two  30  A  cubic 
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structures  (degree  of  polymerization  X  =151).  Thus  we  obtain  the  following 
ranges  for  the  Hildebrand  solubility  parameters  5  (limits  are  standard 
deviations). 

6  =  19.4  ±  1.1  (J/cm3)1/2  113  structures.  X  =  35)  (4a) 

6  =  19.7  ±  0.5  (J/cm3)1/2  (2  structures.  X  =  151)  (4b) 

The  reported  experimental  values  for  the  solubility  parameter  for  PC  are  20.3 
and  20.1  (J/cm3)1/2.3®’31  Hence,  the  theoretical  predictions  provide  quite 
reasonable  results. 

4.2  Conformation  of  the  Individual  Chains  in  Bulk 

There  are  several  aspects  of  the  individual  chain  behavior  in  the  bulk 
that  can  be  investigated  with  the  model  and  compared  with  experiment.  The 
first  of  these  is  the  random  coil  behavior.  As  previously  stated,  in  the  initial 
guess  generation  the  chain  conformations  demonstrate  random  coil  statistics. 
As  with  PP  and  PVC.  the  <r2>0/M  of  the  overall  chain  conformation  in  the 

microstructures  does  not  dramatically  change  during  minimization.17*18  For 
the  fully  minimized  structures,  the  average  <r  2  >0/M  for  the  X  *  35 

microstructures  is  (standard  deviation  in  parentheses)  0.96  ±0.40 
(A  mol  g'1)  sad  for  the  X  >  151  microstructures  the  average  was  0.90  ±  0.14 
(A  mol  g*1).  which  again  compare  well  with  the  reported  experimental  range 
of  0.75  to  1.28.7*8*27*28 

Another  important  area  to  investigate  is  the  influence  of 
intermoiecular  packing  on  the  chain  conformation  with  respect  to  the 
intramolecular  ground  state  conformations  of  the  isolated  molecule.  NMR 
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studies  show  that  site  inhomogeneity  due  to  intermolecular  packing  must  exist, 
while  at  the  same  time  intramolecular  forces  dominate  the  preferred  bulk 
chain  conformations.*  1 5  These  at  first  seemingly  contradictory  Findings  are 
supported  in  our  model.  Figures  3,  4,  and  5  reveal  the  effects  of  intermolecular 
packing  on  the  conformations  of  the  individual  torsion  angles.  In  the  PC 
repeat  unit  (see  Figure  1),  the  conformations  can  be  grouped  into  3  types: 

1.  For  the  two  torsion  angles  and  <p  2  which  determine  the  rotation  of 

the  phenylene  rings  with  respect  to  the  the  isopropylidene  group,  the 
intramolecular  ground  state,  determined  through  calculations  done  on  2.2- 
diphenylpropane,  occurs  when  the  phenylene  rings  are  rotated  out  the  the 
plane  in  a  propeller'  like  conformation.  The  results  discussed  in  (I)  show  that 
the  rotation  of  the  two  angles  is  highly  interdependent.^  3 

2.  The  two  torsion  angles  93  and  95  indicate  the  conformation  of  the 

phenylene  rings  with  respect  to  the  carbonate  group  that  they  bracket.  Using 

diphenylcarbonate  as  a  model  compound  shows  that  the  intramolecular 

ground  states  for  these  torsion  angles  are  not  highly  interdependent  and 
occur  when  the  phenylene  rings  are  rotated  out  of  the  plane  of  the  carbonate 
group  by  about  45°.^  3 

3.  Finally,  94  and  95  determine  the  conformation  of  the  carbonate 

group  locally.  Using  diphenylcarbonate  aa  a  model  reveals  that  the  ground 

intramolecular  energy  states  exist  when  the  carbonate  group  is  in  its 

trans, irons  state,  with  additional  minima  in  the  cis.trans  or  trans.cis  states. ^  3 
The  rotation  of  these  two  angles  exhibit  s  strong  interdependency. 


During  the  analysis,  no  visible  difference  could  be  seen  between  the 
conformations  of  the  18  A  cubes  vs.  the  30  A  cube,  so  the  compilation  of  the 
individual  torsion  angle  conformations  of  ait  15  microstructures  together  is 
compared  with  the  intramolecular  ground  states  and  RIS  states  in  Figures  3,  4 
and  5.  In  Figure  3,  the  background  is  a  contour  map  of  (p  j  vs.  <p 2  computed 

using  2,2-diphenylpropane,  with  x's  located  at  (45°,  45°)  and  at  (135°,  -45°) 
representing  equivalent  ground  state  or  RIS  conformations.  The  energy 
contour  levels  indicate  incremental  increases  of  1  kcal/mole  above  the  ground 
states.  Overlaid  on  the  contour  map  are  the  values  of  the  torsion  angles  that 
exist  in  the  15  microstructures.  The  open  circles  are  the  conformations 
existing  in  the  X  =  35  microstructures,  and  the  filled  circles  are  the 
conformations  existing  in  the  X  =  151  microstructures.  Upon  inspection,  it  is 
obvious  that  the  torsion  angles  are  not  in  their  lowest  energy  conformations 
and  that  there  exists  a  large  variation  of  behaviors.  However,  the  large 
majority  of  the  torsion  angles  occur  in  troughs  of  low  intramolecular  energy 
clustering  around  the  ground  state.  No  conformations  exist  with 
intramolecular  energies  above  about  7  kcal/mole.  This  indicates  that  for  the 
torsion  angles  in  the  isopropylidene  group,  the  intennolecular  packing 
influences  create  a  large  variation  of  torsion  angle  combinations,  but  it  is  the 
intramolecular  energy  contributions  that  determine  the  reference  contour 
giving  which  combinations  can  or  cannot  occur,  in  broad  agreement  with 
NMR  studies.  10-13 

For  torsion  angles  93  and  <p$,  these  two  angles  are  not  closely 
interdependent,  so  the  analysis  can  be  presented  as  a  bar  graph.  In  Figure  4, 
the  individual  torsion  angle  conformations  occurring  in  the  15 
microstructures  are  represented  as  'Dirac  pulses'.  The  intramolecular  ground 
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states  or  RIS  confCTTnaticns  are  located  at  about  -*-45°  and  -45°  out  of  the  plane 
of  the  carbonate  group.  The  figure  shows  that  again  the  conformations 
occurring  in  the  bulk  are  varied,  but  are  none  the  less  centered  around  the 
intramolecular  ground  states.  This  is  in  agreement  with  the  finding  of 
Henrichs  and  Nicely  that  the  phenylene  rings  are  rotated  out  of  the  plane  of 
the  carbonate  group  in  the  bulk.  13  The  intramolecular  energy  barriers  for 
these  torsion  angles  are  much  less  steep  than  the  contours  that  occur  in  the 
isopropylidene  group,  and  therefore  show  more  scatter. 

This  large  scatter  of  torsion  angle  values  becomes  even  more  apparent 
when  studying  the  conformations  of  the  carbonate  group,  94  and  9 5.  In 
Figure  5  the  background  is  a  contour  map  of  torsion  angles  94  vs.  95 
calculated  using  2,2-diphenylpropane.  The  energy  contour  levels  indicate 
incremental  increases  of  1  kcal/mole  above  the  trans, trans  conformation.  The 
ground  states  indicated  with  x's  are  located  at  (0°,  0°)  for  the  trans, traits 
conformation  and  at  (0°.  180°),  (0°,  -180°).  (-180°,  0°)  and  (180°,  0°)  for  the 
trans, cis  and  cis, trans  conformations.  The  carbonate  group  is  much  ’softer’ 
than  the  isopropylidene  group,  resulting  in  larger  areas  of  low  energy  and 
less  steep  energy  well  walls.  Overlaid  on  the  contour  map  are  the  carbonate 
conformations  existing  in  the  15  microstructures.  Again,  the  open  circles  are 
the  conformatioat  existing  in  the  X  *  35  microstructures,  and  the  filled  circles 
are  the  confonattkma  existing  in  the  X  3  151  microstructures.  Immediately  we 
see  that  the  large  majority  of  the  carbonate  groups  are  not  planar  and  that 
there  is  a  wide  scatter  of  conformations,  including  several  with  quite  high 
intramolecular  energies.  The  presence  of  the  high  intramolecular  energy 
conformations  of  the  carbonate  group  indicates  that  the  intermolecular 
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packing  influences  the  conformation  of  the  carbonate  group  more  than  for 
the  phenylene  rings,  which  are  more  dominated  by  intramolecular  forces. 

4.3  Structure  of  the  Glassy  Polymeric  Bulk 

Structural  detail  of  the  intermolecular  packing  can  be  gained  using  the 
bond  direction  correlation  function’  32  which  is  used  to  measure  the 
orientation  correlation  between  bonds.  The  function  is  described  by  the 
equation: 

Sjj  =  (  3<  cos^eij  >  -  l]/2  (5) 

where  9ij  is  the  angle  along  the  chain  between  two  segments  (bond  chords)  of 
a  chain.  The  intermolecular  correlation  function,  SiQter(r).  is  obtained  by 
averaging  over  segments  of  different  images  of  the  parent  chain  with  a  given 
distance  r  between  chord  midpoints.  It  indicates  the  frequency  of  various 
intermolecular  correlations  at  given  distances:  positive  values  indicate  a 
tendency  of  parallel  alignment  of  the  segments,  while  negative  values  reflect 
a  tendency  towards  perpendicular  alignment.  For  PC,  Sinter(t)  was  calculated 
for  the  carbonate  groups  and  is  shown  in  Figure  6.  It  should  be  noted  that  for 
calculations  based  upon  the  sum  of  pairwise  interactions,  there  are  many 
more  pairs  in  oae  of  the  X  a  151  microstructures  than  in  the  X  a  35 
microstructures,  so  the  calculations  done  on  the  two  larger  systems  are 
comparable  to  calculations  done  on  the  sum  of  the  13  smaller  systems.  Upon 
inspection  of  Figure  6,  one  sees  that  there  is  a  tendency  towards  parallel 
alignment  at  a  distance  of  approximately  3  A.  This  is  the  first  indication  that 


17 


there  is  some  trace  of  weak  intermolecular  order  in  the  PC  microstructures, 
albeit  quite  weak. 

Atomic  pair  distribution  functions  can  further  elaborate  the 
intermolecular  structure.  Pair  correlations  functions,  gij(r),  were  calculated 
for  all  possible  pairs.  This  function  is  defined  by  equating  (4*/V)gjj(r)r2dr  to 
the  probability  of  finding  an  i  atom  and  a  j  atom  at  a  distance  between  r  and 
r  +  dr,  where  V  is  the  volume  of  one  microstructure.  For  PC,  there  exists  9 
different  types  of  atoms23.33  when  the  methyl  group  atoms  are  treated 
explicitly  which  therefore  yield  45  distinct  pairs  of  atoms.  The  total  gij(r)  for 
each  pair  is  the  sum  of  intramolecular  and  intermolecular  contributions, 
where  the  intramolecular  contribution  is  obtained  by  'blowing  up'  the  cube  to 
such  a  degree  to  isolate  the  parent  chain  and  thus  count  only  the  interactions 
belonging  to  the  same  chain.  The  intermolecular  contribution  is  then  the 
difference  between  the  total  and  the  intramolecular  gij(r).  Since  not  all  of  the 
45  plots  can  be  shown  here,  only  those  cases  where  interesting  structural 
information  can  be  seen  are  featured.  The  two  sizes  of  microstructures  are 
treated  separately.  Although  there  are  some  visible  small  differences  based  on 
cube  size,  no  trends  can  be  detected. 

In  the  first  atomic  pair  distribution  function,  hydrogen-hydrogen 
(H-H)  shown  in  Figure  7,33  the  functions  of  both  sizes  of  cubes  level  quickly 
off  to  unity,  showing  that  there  is  no  correlation  of  structure  with  respect  to 
hydrogen  atoms.  The  next  pair  distribution  function  to  be  presented  is  the 
carbonyl  oxygen-hydrogen  (0C«0-H)  pair  and  can  be  seen  in  Figure  8.  As 
expected,  there  exists  a  strong  intramolecular  correlation  due  to  the 
hydrogens  on  the  phenylene  rings  flanking  the  carbonate  group,  but  once 
the  intramolecular  contribution  is  removed,  no  significant  correlation  is 
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found  The  main  chain  oxygen-hydrogen  (O-H)  distribution,  which  is  not 
shown  here,  demonstrates  the  same  amorphous  behavior. 

Some  of  the  pair  distribution  functions  do  show  some  intermoiecuiar 
correlations.  The  isopropylidene  carbon-isopropylidene  carbon  (Ca,-Cal).  the 
carbonate  carbon-carbonate  carbon  (Cc=0-Cc=0),  the  carbonate  carbon- 
carbonate  oxygen  (Cc=0-Oc=0),  and  the  carbonate  carbon-main  chain  oxygen 
(Cc=0-0)  pair  distribution  functions,  featured  respectively  in  Figures  9  to  12, 
all  show  some  intermoiecuiar  correlations  in  both  sizes  of  microstructures. 
These  correlations  are  not  particularly  strong,  but  do  suggest  that  some  areas 
on  the  structure  are  not  strictly  amorphous,  as  was  the  case  with  PP.17  The 
correlations  found  in  the  (Cal-Cal)  pair  distribution  function  shown  in  Figure 
8  indicate  that  there  is  some  weak  intermoiecuiar  preference  for  distances  of 
about  5  to  6  A,  which  could  be  an  indication  of  the  areas  of  'enhanced  order1 
found  using  experimental  wide  angle  scattering  techniques.  ^  *6.9 

Other  structural  features  were  also  examined.  The  spatial  uniformity  of 
the  structures  was  studied,  and  indicated  that  there  are  some  'large' 
inhomogeneities  in  PC  in  the  form  of  empty  space  in  the  box  on  the  order  of  3 
to  3  A  in  diameter.  A  more  rigorous  analysis  of  the  empty  space  was  performed 
by  Arizzi  and  Suter.3*  and  demonstrated  that  the  degree  of  heterogeneity  in  PC 
was  not  a  function  of  the  cube  size.  This  feature  along  with  the  other  analysis 
presented  here  dhow  few  differences  between  the  two  size#  of  micro  structures, 
indicating  that  the  system  is  already  'converged'  at  the  smaller  cube  size,  and 
that  both  sizes  can  be  used  for  future  simulations  of  the  structural  detail  of  PC. 


19 


V.  CONCLUSIONS'* 

The  analysis  of  the  PC  microstructures  show  that  the  PC  chain  exhibits 
random  coil  behavior  in  the  bulk.  Both  the  intramolecular  and  intermolecular 
contributions  determine  the  individual  chain  conformations,  with  the 
intramolecular  forces  dictating  what  conformations  are  possible  and  which 
are  not.  The  intermolecular  forces,  especially  for  the  carbonate  group, 
contribute  to  the  large  variation  of  conformations  that  are  centered  around 
the  intramolecular  ground  states.  On  the  whole,  the  structure  of  PC  is 
amorphous,  although  some  weak  intermolecular  correlations  on  a  very  small 
scale  do  exist.  The  features  studied  show  excellent  agreement  with  experiment, 
indicating  good  quality  of  the  simulated  PC  microstructures  and  their 
usefulness  for  future  simulations. 
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Table  1 


Calculated  Energies  Upar*  AU tails*  and  U  for  PC 


Upar  (kcal/mol  of  structures)4  56.6  ±  16.1  219.3  ±  31.2 

AU  tails  (kcal/mol  of  structures^  22.6  ±  2.6  -121.3  ±6.1 

U  (kcal/mol  of  structures)4  -260.3  ±  20.1  -1146.5  ±  32.7 


4  ^he  range  reported  is  the  standard  deviation. 

b  The  effective  range  for  the  complete  amorphous  and  for  the  liquid  like  order 
approximations  is  reported. 
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REPEAT  UNIT  CHAIN  END 


showing  ihe  virtual  bonds  and  their  associated  torsion  angles 


2.  A  minimized  PC  microstruclure  of  X  ^  35.  In  (a)  the  atomic  radii  have  been  reduced 
for  clarity  while  in  (b)  the  actual  atomic  radii  have  been  included 


3.  The  intramolecular  vs.  the  intermolecular  influence  of  the  pfienylene  ring 
conformations  in  the  isopropyiidene  group.  The  background  is  a  contour  map  of  vs. 

<P2  computed  using  2,2-diphenyipropane,23  with  x's  located  at  (45*.  45*)  and  at 
(135*,  -45’)  representing  the  ground  state?  or  RIS  conformations.  The  energy  contour 
levels  indicate  incremental  increases  of  1  kcai/mole  above  the  ground  states.  Overlaid  on 
the  contour  map  are  the  values  of  the  torsion  angles  occurring  in  the  15 
microstructures.  The  open  circles  are  the  conformations  existing  in  the  X  •  35 
microstructures,  and  the  filled  circles  are  the  conformations  existing  in  the  X  •  151 
microstructures. 
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4  Distribution  ol  torsion  angles  occurring  in  the  15  microstructures  The 
intramolecular  ground  stales  or  RlS  conformations  are  located  at  about  ♦  45u  and  -  45 
and  are  equally  populated  m  the  initial  guess  structures. 
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5.  The  intramolecular  vs.  the  intermodular  influence  for  the  carbonate  group 
conformations.  The  background  is  a  contour  map  of  torsion  angles  94  vs.  95  calculated 

using  diphenyipropane.23  The  ground  states,  located  at  (0*.  0*)  for  the  trans.lrans 
conformation  and  at  (0°,  180*),  (0*.  - 1 80°) ,  (-180*,  0°)  and  (180*.  0*)  for  the 
'rans.Cis  and  cis.trans  conformations,  are  indicated  with  x's.  The  energy  contour  levels 
indicate  incremental  increases  of  1  kcal/moie  above  the  trans.lrans  conformation. 
Overlaid  on  the  contour  map  are  the  carbonate  conformations  that  occur  in  the  15 
microstructures.  The  open  circles  are  the  conformations  existing  in  the  X  •  35 
microstructures,  and  the  filled  circles  are  the  conformations  existing  in  the  X  •  151 
microstructures. 


6.  The  intermolecular  average  bond  direction  correlation  function  Sjnter(0  calculated 
for  the  carbonate  groups. 


Pair  Distribution 
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9.  Isopropylidena  carbon-isopropylidena  carbon  (Ca,-Cal)  atomic  pair  distribution 
function. 


Radial  Distance  (A) 


1 0.  Carbonate  carbon-carbonate  carbon  (Cc*°-Cc*°)  atomic  pair  distribution  function. 


Radial  Distance  (A) 


1 1 .  Carbonate  carbon-carbonate  oxygen  (Cc-0-Oc"°)  atomic  pair  distribution  Junction. 


12.  Carbonate  carbon-main  chain  oxygen  (Cc*°-0)  atomic  pair  distribution  function. 


